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Schwinger-Dyson equation, we estimate the individual non-perturbative contribution of the "one- 
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the gluon propagator is affected by the removal of the ghost loop, and what are the consequences 
on the phenomenon of gluon mass generation. 
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1. Introduction 

In the last few years our understanding of the infrared (IR) behavior of the fundamental QCD 
Green's function has improved substantially. Putting together the information obtained though 
various non-perturbative methods, such as lattice simulations [jlj [| ^ Q], Schwinger-Dyson 
equations (SDEs) [[7], ||, [J, functional methods [ JTOj , |Tl[ ], and algebraic techniques [13, 14], it is by 



now well-established that, in the Landau gauge, the gluon propagator and a ghost dressing function 
are finite in the IR (in d = 3,4). [[7], ||]. Evidently, these results support the gluon mass generation 
picture proposed by Cornwall several years ago JOl), disfavoring the so-called "ghost-dominance" 



picture of QCD |[1^, |17[ |, whose theoretical cornerstone has been the existence of a divergent ("IR- 
enhanced") ghost dressing function. 

However, the finiteness of the ghost sector does not imply necessarily that the ghost contribu- 
tion has been relegated to a marginal role in the QCD dynamics. In fact, compelling evidence to 
the contrary has emerged from detailed studies of the gap equation that controls the breaking of 



chiral symmetry and the dynamical generation of a constituent quark mass Ql 8| , |19| ] . Specifically, a 
detailed study has revealed that the proper inclusion of the ghost sector into the quark SDE is cru- 
cial for obtaining quark masses of the order of 300 MeV in the presence of finite gluon propagator 



[19] 



Given the importance of the ghost sector for the dynamical generation of a constituent quark 
mass, the main purpose of this talk is to ask whether a similar situation applies in the case of the 



dynamical generation of an effective gluon mass [20]. 



2. Disentangling the "one-loop dressed" ghost contributions 

In what follows we will work within the specific framework provided by the synthesis of the 



pinch technique (PT) [[15|, |21|, |SJ, |22j, |2J] with the background field method (BFM) [|24 

We start by recalling that within the PT-BFM scheme the gluon self-energy, Yl llv {q), is giving 
by the sum of the diagrams represented by Fig. [IJ, i.e. 

10 

n^) = I>r. (2.i) 

1=1 

Exploiting the well-known blockwise transversality properties of the PT-BFM YV iV {q), it is 
possible to separate the transversal contribution of the one-loop dressed ghost diagrams, repre- 
sented by the diagrams (a?,) and (04), and to be denoted by 

nr(q) = (a 3 r + (a 4 y v ; (2.2) 

where (#3) and (a 4 ) are given by 

O3W = -g 2 C A Jr^\k,q,-k-q)D(k)D(k + q)f v (k + q,-q,-k), 

(a 4 ) MV = 2g 2 C Agflv jD(k). (2.3) 
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Figure 1: The SDE corresponding to the PT-BFM gluon self-energy Y\ a ^ v {q). The graphs inside each box 
furnish an individually transverse contribution. White (black) circles denote full propagators (vertices). 



In the equations above, D ab (q 2 ) = 8 ah D{q 2 ) denotes the full ghost propagator, defined in terms of 
the ghost dressing function F(q 2 ) as 

D( q 2 ) = ^, (2.4) 

while represents the three-particle vertex describing the interaction of the background gluon 
with a ghost and an antighost, with (all momenta entering) 

iT c b Al 1 A r ^^P) = Sf acb ^ii{ r ^^PY rj, 0) (r,<7,p) = {r-p)^. (2.5) 

Finally, Ca is the Casimir eigenvalue of the adjoint representation [Ca = N for SU(N)], and we 
have introduced the <i-dimensional integral measure (in dimensional regularization) according to 



d a k, (2.6) 



" (2n) d . 

with jU the 't Hooft mass, and e = 4 - d. Then, by virtue of the PT-BFM Ward identity 

iq^f^r,q,p)=D- l (r)-D-\p), (2.7) 
it is immediate to establish the transversality of Tic v (q), namely p2[ ] 

q^ v (q)=0. (2.8) 
It is convenient for our purposes to decompose the full self-energy Yl llv (q) as 

iP v ( 9 )=nr(?)+nr(?), (2-9) 

where IT^ v (q) denotes the sum of the remaining subsets of diagrams in Fig. |[ i.e., both the gluon 
one- and two-loop dressed diagrams, as well as two-loop dressed ghost diagrams, 

10 

nr(*)= L(«<r- (2.io) 

i=\ 
#3,4 
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Notice that due to the special Ward identities satisfied by the PT-BFM vertices, IT^ (q) is also 
transverse [0, 23]. 



Using Eq. ( |2.9[ ), the SDE for the full gluon propagator in the Landau gauge of the PT-BFM 
scheme assumes then the form [Q] 



A-V)^ V (<7) 



q 2 P^(q)+i[nr(q)+nr(q)] 



[l+G(q 2 )Y 



(2.11) 



where the gluon propagator A^ v (q) is defined as 



A fiv(q) = -iA(q 2 )Pnv(q); P^iv{q) = 



q^q v 



(2.12) 



The function G appearing in ( |2.11| ) is the form factor associated with g^ v in the Lorentz decompo- 
sition of the auxiliary two-point function A, given by [[7], |25|] 



A Mv(^) = ~ig 2 C A jA°(k)D(q - k)H va (-q,q - k,k) 



^ V G(, 2 ) + ^L(, 2 ) 



(2.13) 



Notice that the auxiliary function H is related to the (conventional) gluon-ghost vertex by the 
identity 

r fi( r ,<l,P) = -p V H v ^ l (p,r 1 q), (2.14) 



and that, in the (background) Landau gauge, the following all order relation holds [|2g, [270 

F-\q 2 ) = \ + G{q 2 )+L{q 2 ). (2.15) 
Now, let us return to Eq. ( [2.1 1[ ), and define in a completely analogous way the quantity A,-(q 2 ), 

q 2 P^{q) + iU^{q) 



given by 



A-\q 2 )P^(q) 



[l+G(q 2 )Y 



(2.16) 



Evidently, A r represents the propagator obtained by subtracting out from the full propagator A the 
one-loop dressed ghost contributions. Then, taking the trace of both Eqs. (2.1 1) and (2.16), defining 
the trace of lie v (q) as 

n c ( q 2 ) = n^( q ), (2.17) 



and solving for A r , we arrive at 

A r {q 2 )=A(q 2 )\\ 



iA(q 2 )Yl c ( 



-1 



{d-\)[l+G(q 2 



(2.18) 



which represents our master formula. 

In order to obtain the behavior of the propagator A r {q 2 ) from Eq. (2.18) we will in the next 
sections (i) identify the full gluon propagator A(q 2 ) with that obtained from the lattice, and (ii) 
determine nonperturbatively the quantity II C from Eqs. ( |2T3| ) and ( |2.17| ), and evaluate it numerically 
using as input the lattice results for the ghost dressing function F{q 2 ). 
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3. The nonperturbative calculation of Ti c (q 2 ) 

The first step in the calculation of the non-perturbative quantity Yl c (q ) is the definition of an 
Ansatz for the fully-dressed ghost vertex , appearing in graph (as ) of Eq. ( |2.3| ) which satisfies 



the crucial Ward identity of Eq. ( |2.7| ). This task can be accomplished with the help of the "gauge- 



technique" [Eq] which reconstruct the vertex by "solving" its Ward identity. Following the same 



steps of the derivation presented in Q29Q for the scalar QED vertex, the Ansatz for the fully-dressed 
ghost vertex , reads, 

f^p) = Act [£>~V)-£>~V)] , (3.i) 



which evidently satisfies Eq. ( [277| ) when contracted with . Obviously the "gauge technique" 
leaves the transverse (automatically conserved) part of the vertex undetermined, which, on general 
grounds, is expected to be subleading in the IR [28, 3CQ . 



Substituting ( |3.1| ) in the first equation of (2.3) and taking the trace, it is relatively straightfor- 
ward to obtain the result 

U c (q 2 )=g 2 C A [AT{q)-q 2 R{q)} , (3.2) 

Run - i D ( k +<i)- D ( k ) 



where 



k 



m= l ^tT +iim. (33) 



k 



(k + q) 2 -k 2 ' 
2 D(k + q)-D(k 
Ik (k + q) 2 -k 2 ' 2 

To further evaluate Yl c (q 2 ), we must invoke the so-called "seag ull-identity" [fl]], 

valid in dimensional regularization, which enforces the cancellations of all seagull-type of diver- 
gences. Notice that, in the limit q — > 0, the term q 2 R(q) vanishes, and so does T(q), since 

where in the last step we have employed Eq. (|3~/i|), with f(k 2 ) — > D(k 2 ). Employing this result, 
follows immediately from Eq. (3^) that n c (0) = 0. Using the fact that A _1 (0) = m 2 (0), we con- 



clude that the one-loop dressed ghost diagrams (03) and (a^) do not contribute directly to the value 
of dynamical gluon mass at zero momentum transfer. The easiest way to appreciate this is by re- 
calling that the mechanism responsible for endowing the gluon with a dynamical mass relies on 
the presence of massless poles in the nonperturbative tree-gluon [the black circle in graph (ai) of 
Fig. whereas the ghost vertex has the usual structure [note the absence of poles in the Ansatz of 



Eq. (£TD] V£3- 

In addition, notice that when d = 4, R(q) is ultraviolet divergent, and must be properly renor- 
malized, by introducing the appropriate wave-function renormalization constant. 



5 



Impact of ghost loops on dynamical gluon mass generation 



Arlene C. Aguilar 



The (subtractive) renormalization must be carried out at the level of (2.11). Specifically (set 
ting directly d = 4), 

Z A q 2 + Un r (q)+U c (q)] 



i-l/ 



(3-6) 



[i+G(<7 2 )r 

where the renormalization constant Za is fixed in the MOM scheme through the condition A -1 (ju ) = 
Applying this condition at the level of Eq. ( |3.6[ ) together with Eq. ( |2.15| ), and using the fact that the 

nail 
F 

Z A 



function L(x) is considerably smaller than G(x) in the entire range of momenta, (so that we can use 



the approximation 1 + G(/i 



-V) 
1 - 



1) allows one to express Za as 



3/1 



2[n,-Gu)+n c Gu)] 



(3.7) 



Finally, substituting Eq. ( |3.7[ ) into Eq. Q3.6| ), and defining (in a natural way) the renormalized 

A~ l (q 2 ) as 



A _i ( 2) = ^ 2 + i[n r (g)-(g 2 /M 2 )n,(M)] 



[l + G(q 2 )] 2 



the renormalized version of the master formula (J2. 1 8p will read 

A r ~V) 



k-li 



i [n c ( 9 )-( 9 2 /M 2 )n c (M)] 



(3.8) 



(3.9) 



Evidently ( |3~9^ ) is obtained from ( 2. 1 8 ) by replacing A 1 {q 2 
and H c (q) — > n c ft(^), where 

n c , s («?)=n c (^)-(^ 2 / i u 2 )n c (At). 



1 + G(</ 2 )] 2 

A^ 1 ^ 2 ) ("/?" for "renormalized"), 



(3.10) 



For the ensuing numerical treatment of R(q) and T{q) carried out in the next section, it is 
advantageous to have the crucial property T(0) = a priori built in, in order to avoid possible 
deviations due to minor numerical instabilities. To that end, we introduce the quantity T 

~D(k + q)-D(k) dD(k)~ 



T(q)=T(q)-T(0) 



(k + q) 2 -k 2 dk 2 



(3.11) 



which has the property of ensuring (by construction) that T (0) = 0, while, at the same time, coin- 
ciding with the original T for all momenta q. 

In addition, it is convenient to re-express R(q) and T(q) in terms of the ghost dressing function. 
Using Eq. (2.4), after some elementary algebra, one obtains 



R(q) 



F(k) r F(k + q)-F(k) 

p(k + q) 2+ J k k 2 [(k + q) 2 -k 2 y 



T{q) 



F{k + q)-F(k) dF{k) 
{k + q) 2 -k 2 ~W 



;(3.12) 



note that the angular integration of the first term in R can be carried out analytically for any value 
of the space-time dimension d. 

Finally, note that up until this point we have been working in Minkowski space. To make the 
transition to Euclidean space, we must employ the usual rules. Specifically, we set J k = i J k and 
q\ = -q 2 , and use that A E (<? 2 ) = -A(-<7 2 );F E (g 2 ) = F(-^r 2 );G E (^r 2 ) = G(-ql), suppressing the 
subscript "E" in what follows. 
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Figure 2: Left panel: Lattice result for the ST/ (3) gluon propagator, A(q), in d = 4, renormalized at 
ji — 4.3 GeV. The continuous line represents the fit given by Eq. (4.1). Right panel: The SU(3) ghost 
dressing function, F(q 2 ), renormalized at the same point, fi = 4.3 GeV; the solid line corresponds to the fit 
given by Eq. (4.2). 



4. Numerical Results 



We will now proceed to perform the numerical analysis. Using the available lattice data on the 



ghost dressing function F, we evaluate the terms R and T given in Eq. (3.12), and combine them 
following the Eqs. (3.2) and ( 3.10| ) to obtain the (renormalized) ghost contribution to the gluon 



self-energy Tl c . Finally, we construct A,- using (2.18) and the lattice results available for the gluon 
propagator A. This exercise is carried out for two different cases: d = 4, N = 3, and d = 3, N = 2. 

In Fig. H we show the lattice results for the four-dimensional SU (3) gluon propagator A(q 2 ) 
(left panel), and the corresponding ghost dressing function F(q 2 ) (right panel), obtained from [g], 
and renormalized at = 4.3 GeV. 



As has been discussed in detail in the literature |32j, |33|], both sets of data can be accurately 
fitted in terms of IR-finite quantities. More specifically, for the case of A(q 2 ), we have proposed a 
fit of the form [j33h 



A-V)=MV) + <7 : 



■In 



l3CAgj ,_fq 2 + Pl M 2 (q^ 
967T 2 



M\q 2 ) 



q 2 + p 2 m 2 



(4-1) 



Notice that in the above expression, the finiteness of A~ l (q 2 ) is assured by the presence of 
the function M 2 (q 2 ), which forces the value of A -1 (0) = M 2 (0) = m 2 ,/ p 2 . The continuous line on 
the left panel of Fig. ||| corresponds our best fit, which can be reproduced setting mo = 520 MeV, 
g\ = 5.68, pi = 8.55 and p 2 = 1.91. 

The SU{3) lattice data for F(q 2 ), shown in the right panel of Fig. ||, will be fitted by the 
following expression 



F-\q 2 



1 + 



9C A g 2 1 Jq 2 +p 3 M 2 (q 2 ' 



4 487T 2 



In 



M 2 (q 2 ) 



q 2 + p 4 m 2 ' 



(4.2) 



with the parameters given by mo = 520 MeV, g 2 = 8.65, P3 = 0.25 and p 4 = 0.64. Notice that the 
M(q 2 ) has the same power-law running as the one reported in Refs [34, 35, [36|]. 
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Figure 3: Left panel: Numerical evaluation of the ghost contribution Tl c (q) to the gluon propagator using 
as input our best fit for the d = 4, N — 3 ghost dressing lattice data. Right panel: The removal of the one- 
loop dressed ghost contribution from the (lattice) gluon propagator results in a diminished "swelling" in the 
momentum region below 1 GeV 2 . 



The only missing ingredient for the actual nonperturbative determination of II C , and therefore 
A r , is the value of a s = g 2 j4%. Instead of choosing a single value for a s , we will use the physi- 
cally motivated range of values [0.20,0.29], which will furnish a more representative picture of the 
numerical impact of the ghost corrections on the gluon propagator. 

The results obtained for the renormalized R and T, after substituting into the corresponding 
formulas our best fit for F, given by Eq. are shown on the left panel of Fig. ||, together with 
the combination q 2 R — 4T, which appears on the rhs of Eq. It is clear that the contribution 

of the term 4T is rather negligible; in a way this is to be expected, given that this term vanishes 
identically in perturbation theory (for all values of q), and vanishes nonperturbatively at the origin. 



Next, we use these results to construct IT C , given in Eq. (3.1C), and finally A r , expressed by 



Eq. ( |3T9| ) (Fig. g right panel), using both values of a s , namely a s = 0.29 (blue dashed line) and 
a s = 0.20 (magenta dashed-dotted line). 

We then see that the net effect of removing the ghost contribution is to suppress significantly 
the support of the gluon propagator in the region below 1 GeV 2 . Higher values of a s increase the 
impact of the ghost contributions, but only slightly, as can be seen on the right panel of Fig. ||. As 
we will see in the next section, this "deflating" of the gluon propagator in the intermediate region 
of momenta, produced by the removal of the ghost contributions, has far-reaching consequences 
on the generation of a dynamical gluon mass. 

Now we will repeat the same exercise using the lattice results for d = 3 and N = 2. Let us start, 
as in the previous case, by showing in Fig. ^ the lattice results [Q] for the three-dimensional gluon 
propagator A(q) (left panel) and the ghost dressing function F(q) (right panel). Both A(q) and 
F (q) saturate in the deep IR region, and can therefore be fitted by means of IR finite expressions. 

In the case of the gluon propagator, an accurate fit is giving by 

A( q )=A^[-( q - qo f/ W ] + a+b 1 q + cq2 , (4.3) 
where the fitting parameters are A = 0.49, q = 0.11, w = 0.37, a = 0.43, b = -0.85, and c = 1.143. 
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Figure 4: Left panel: Lattice results for the SU (2) gluon propagator in d = 3. The continuous line represents 
our best fit to the data obtained from Eq. (4.3). Right panel: Lattice data for the SU(2) ghost dressing function 
F(q) in 3 dimensions; the solid line corresponds to the best fit given by Eq. (4.4). 



For the ghost dressing function, we use the following piecewise interpolator 

= —T-> 2 > for I 2 ^ 3 and = 1 + — 7-2 ' for f > 3 ' (4 " 4) 

a + bq + cq L eq + q L 

with fitting coefficients a = 0.19, b = 0.61, c = —0.14, d = 0.63 and e = 0.26 obtained by requiring 
the function to be continuous at q 2 = 3. 



Next, substituting the results presented on the left panel of Fig. [5| into Eqs. ( p.2[ ) and ( |2.18| ), 
and using g = 1 .208, we compute IT C and A r . On the right panel of Fig. |5l we compare the residual 
propagator A r (blue dashed line) with the full propagator A(q). Clearly, the effect in the tridimen- 
sional case is even more pronounced: the ghost contribution completely dominates over the rest, 
determining to a large extent the overall shape and structure of the propagator. 




Figure 5: Left panel: Numerical evaluation of the ghost contribution n c to the gluon propagator using as 
input our best fit for the d = 3, = 2 ghost dressing lattice data. Right panel: The result of removing the 
one-loop dressed ghost contribution from the gluon propagator in d = 3. The effect is much more dramatic 
than in the d = 4 case, since all the structure is determined by the ghost contribution, while A r has the sole 
(but crucial!) role of rendering the propagator finite at q = 0. 
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□ - 



0.0 0.4 0.8 12 1.6 2.0 o.o 0.5 1.0 1 5 2.0 2 5 3 3.5 4.0 

q [arbitrary units] q [arbitrary units] 

Figure 6: The kernel Jt^/j^ of Eq. (5.2) constructed out of the lattice propagator A (left panels) and the 
ghost-less propagator A,, (right panels) for the d = 4 N = 3 (top row), and d = 3 N = 2 (bottom row) cases. 
The insets show in each case the shape of the propagator used to evaluate the kernels. 



5. The effects of the ghost loop in the mass equation 

In the previous section we have studied how the subtraction of the ghost contributions affects 
the profile of the gluon propagator. However, as we will see now, the effects goes way beyond a 
simple change in the overall propagator shape, modifying its salient qualitative characteristics, and 
in particular the generation of a dynamical gluon mass. 

To establish this, we use the q — > limit of the equation describing the behavior of the dynam- 
ical mass equation, i.e. 

m2 (°) = r^ 2 (ymMy), (5.i) 

d(4w)sr(|) l + G(0) Jo 

with the kernel J(fd;N given by 

je d . N (y)=yi- l A(y)\yA(y)}'. (5.2) 



Since the constant multiplying the integral is positive, the negative sign in front of Eq. ( p.l| ) 
tells us that the required physical constraint m 2 (0) > can be fulfilled if and only if the integral ker- 
nel J%i;N (constructed solely out of the gluon propagator) displays a sufficiently deep and extended 
negative region at intermediate momenta [p^]. 

In the left panels of Fig. |6] we plot the kernels J% ; jv obtained from the lattice data for the cases 
d = 4, N = 3 (top row), and d = 3 N = 2 (bottom row); both display the characteristic negative 



region that allows, at least in principle, the existence of solutions of Eq. (5.1), furnishing a positive 
value for m 2 (0). 
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On the other hand, the situation changes substantially once the ghost loop is removed, in which 
case the kernels <%d;N must be constructed from A r (right panels of the same figure). For d = 4 one 
observes a shift towards higher qs of the zero crossing, and a correspondingly suppressed negative 
region; even though this is not sufficient to exclude per se the existence of a physical solution to the 
mass equation (^jj), a thorough study of the approximate equation derived in [ |32| ] reveals that no 
physical solution may be found. The d = 3 situation is even more obvious: the highly suppressed 
negative region present in this case cannot support solutions of (5T) with ra 2 (0) > 0, thus leaving 
as the only possibility the trivial m 2 = solution. 

The main conclusion one can draw, therefore, is that the ghosts play a fundamental role in the 
mechanism of dynamical gluon mass generation, since the failure to properly include them results 
in the inability of the theory to generate dynamically a mass for the gluon. 

6. Conclusions 

In this talk we have presented a study of the impact of the ghost sector on the overall form of 
the gluon propagator in a pure Yang-Mills theory, for different space-time dimensions (d = 3, 4) 
and SU (N) gauge groups (N = 2, 3). 

The suppression of the gluon propagator induced by the removal of the ghost-loops has far- 
reaching consequences on the mechanism that endows gluons with a dynamical mass, associated 
with the observed IR-finiteness of the gluon propagator and the ghost-dressing function. Specif- 
ically, using a recently derived integral equation controlling the dynamics of the (momentum- 
dependent) gluon mass, we have demonstrated that when the reduced gluon propagators are used 
as inputs, the corresponding kernels are modified in such a way that no physical solutions may be 
found, thus failing to generate a mass gap for the pure Yang -Mills theory pOj). Instead, as has been 



shown in [32], the use of the full gluon propagator in the same equation generates a physically 



acceptable gluon mass. 
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